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ABSTRACT 

We present Very Large Array (VLA) high resolution observations of the NH3(1,1) 
and NH3(2,2) molecular transitions towards the high mass star forming region AFGL 
437. Our aim was to investigate if the poorly collimated CO molecular outflow pre- 
viously detected in the region is the result of a projection effect, with no intrinsic 
bipolarity, as suggested by Gomez et al. We complemented our observations with ra- 
dio continuum archived data from the VLA at 2 and 3.6 cm, and with unpublished 
public data at 450 /xm taken with Submillimetre Common-User Bolometer Array at 
the James Clerk Maxwell Telescope. Ammonia emission was found mainly in three 
clumps located at the south and east of the position of the compact infrared cluster 
of AFGL 437, where the CO outflow seemed to have its origin. One of the NH3(1,1) 
clumps coincides with the maximum of NH3(2,2) and with a local peak of emission at 
450 /im. A near infrared source (sll) is also found at that position. Our continuum 
map at 2 cm shows extended elongated emission associated with the infrared source 
AFGL 437W. This elongated morphology and its spectral index between 3.6 and 2 
cm (~ 0.4) suggest the presence of a jet in AFGL 437W. We suggest that several 
molecular bipolar outflows may exist in the region. The observed CO outflow would 
be the superposition of those individual outflows, which would explain its low degree 
of collimation observed at larger scales. 
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1 INTRODUCTION 

Mass-loss phenomena are one of the best-known manifesta- 
tions of the star formation process. In our current under- 
standing of the formation of low-mass stars via accretion 
(e.g., McKee & Ostriker 2007), highly collimated jets are 
necessarily present in the first stages of protostellar evolu- 
tion, coeval with the formation of circumstellar disks. Jets 
could be the agent releasing angular momentum excess, so 
that the protostar can continue accreting material from its 
environment. 

Evidence for highly collimated mass-loss is widespread 
in low-mass young stellar objects (YSOs): e.g., jets traced by 
Herbig-Haro objects, radio continuum emission, or masers. 
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as well as bipolar molecular outflows. However the case for 
high-mass stars (M J5 8 M0) is less clear. Energetic mass 
loss is indeed present in those sources, althoug h it is, in 
gener al, less collimated than in low-mass objects (|Wu et al.l 
|2004 ). Highly collimated jets seems to be restricted to the 
earliest phases of the evolution (< 10^ years; Shepherd 2005) 
of massive YSOs. 

This is important from a more global perspective, since 
it is not yet clear whether high-mass stars form via accretion, 
like their low-mass counterparts (e.g., Yorke & Sonnhalter 
2002; McKee & Tan 2003; Krumholz et al. 2009), or by co- 
alescence of lower-mass objects (Bonnell, Bate & Zinnecker 
1998). The presence of circumstellar disks and collimated 
outflows are key in gredients of the a ccretion scenario. The 
detection of disks ( Patel et al.l l2005h and jets (e.g., Marti, 
Rodriguez & Reipurth 1993; Rodriguez et al. 1994; Davis 
et al. 2004; Patel et al. 2005) in several high-mass objects 
clearly indicate that formation via accretion is possible in 
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Figure 1. Integrated intensity map of the main hyperfine component of the NH3(1,1) (left) and NH3(2,2) (right) Unes. The lowest 
contour levels and the increment step are 3 times 8.7 mjy km , the rms of the map (beam= 8", shown at the bottom left-hand 
corner o f each map). Ammonia clumps are lab elled as CI, C2, and C3. Crosses mark the position of the infrared sources detected with 
Spitzer llKumar Dewangan fc Anandaradl2010l) . of which the better studied ones are labelled. Axes are right ascension and declination 
in equinox J2000. 



these objects, but we still do not know if this process is of 
general application. 

Here we will study the AFGL 437 region, which hosts an 
interesting molecular outflow. The AFGL 437 region com- 
prise s a cluster of at least ^ 20 YSOs JWeintraub fc Kastnerl 
1 19961 : E umar Dewangan fc Anandarad [2010l ). although in- 
frared images are dominated by a compact (~ 15") central 
cluster of four sources (named AFGL 437N, S, E, and W 
by Wynn- Williams et al. 1981), which seem to have recently 
emerged from the near side of the molecular cloud. AFGL 
437 shows clear signs of ongoing star formation, such as 
the presence of radio continuum emission (probably trac- 
ing ultracompact HII regions associated with sources W 
and S, Wynn-Williams et al. 1981; Torrelles et al. 1992; 
Kurtz, Churchwell fc Wood 1994), water masers (close to 
sources N and W, To rrelles et al. 1992), and a molecular 
outflow traced by CO (jGomez et al ]Ll992). At least sources 
W and S are thought to be massive stars, of ear ly B type 
IWvnn- Williams et al. I [T98ll : I Torrelles et al.lll992l ). 

The molecular outflow, observed w ith single dish at 13" 
angular resolution (|G6mez et al.|[l992i ). is roughly oriented 
in the north-south direction, but it shows a very low degree 
of coUimation, with high-velocity CO emission completely 
surrounding the central cluster. Interestingly, all observed 
CO isotopes show the same distribution of blue- and red- 
shifted gas. Given the location of th e YSO s close to the edge 
of the parental cloud, IComez et al.l (| 19921 ) suggested a pos- 
sible interpretation for the nature of the molecular outflow, 
alternative to more classical models: if the winds from the 
YSOs, that can be isotropic at the origin, shock obliquely 
against the walls of the cavity opened in the cloud, they can 
produce a laminar flow of molecular gas along these walls. 
When this flow is observed under a particular angle of view, 
it may give the appearance of a bipolar molecular outflow. 
The observed bipolarity would then be a projection effect, 
not an intrinsic characteristic of the mass-loss itself. An al- 



ternative explanation for the morphology of the molecular 
outflow would be the superposition of several outflows ex- 
cited by different sources in the cluster. 

The IR observations sh owed two different s ources 
towards AFGL 437N (jRavner fc McLeanI Il987l : 
IWeintraub fc Kastneil Il996l ). one of which (WK 34) is 
associated with a bipolar polarized nebula (Weintraub 
& Kastner 1996; Meakin, Hines fc Thompson 2005) 
oriented N-S (i.e., roughly in the same direction as the 
CO outflow) w ith a centrosymmetri c pola risation pat- 
tern. This led IWeintraub fc KastneJ (| 19961 ) to suggest 
that WK 34 is the main driving source of the molecu- 
lar outflow. How ever, the na t ure o f this source is still 
uncertain: while , Meakin et al.l (|2005l ) argued that it is a 
low- m ass YSO, recent SEP modelling inc luding Spitzer 
data iKumar Dewangan fc Anandarad l2010i ) suggest that 
it could be massive, but of young age, with effective 
temperature still not sufficient to create an HII region. On 
the other hand, high-resolution infrared observations show 
elongated ( monopolar) emission associated with source 
AFGL 437S (| Alvarez et al.ll2004l ). This suggests that more 
than one source in the region could be undergoing mass 
loss. 

In this paper we present high resolution ammonia ob- 
servations of the AFGL 437 region, designed to test whether 
the bipolarity of the outflow is as a projection effect (as sug- 
gested by Gomez et al. 1992) or if, alternatively, the observed 
outflow is the superposition of individual outflows from dif- 
ferent sources. In the case of a projection effect, we should 
expect the ammonia to show a spatio-velocity pattern simi- 
lar to that of CO isotopes, with blueshifted gas to the south 
and redshifted one to the north. On the other hand, if this 
is a more classical example of molecular outflow, we expect 
ammonia emission to peak near the exciting source of the 
outflow, and possibly tracing a disk/toroid elongated in the 
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E-W direction, i.e., perpendicular to the collimation axis 
(e.g., Estalella et al. 1993; Wiseman et al. 2001). 

We complement the ammonia observations with 
archival data of continuum emission at centimetre and sub- 
millimetre wavelengths to further study the nature of the 
sources in this region. 



2 OBSERVATIONS 

We observed the (1,1) and (2,2) inversion transitions of the 
NH3 molecule (rest frequencies 23694.496 and 23722.634 
MHz, respectively) towards the AFGL 437 region, using 
the Very Large Array (VLA) of the National Radio As- 
tronomy Observatorjo in its D configuration on 2004 Au- 
gust 16 (project AG665). The phase centre was located at 
R.A.(J2000) = 03'^07'"23.7=, Dec(J2000) = 58°30'50". We 
observed both lines simultaneously in dual circular polar- 
isation, using the 4IF mode of the VLA. This allowed us 
to sample each line with 63 spectral channels of 48.8 kHz 
(0.62 km s"^) width, centred at Vlsr = -39.4 km s"^ The 
flux calibrator was J0137+331, with an assumed flux den- 
sity of 1.05 Jy using the latest VLA values (1999.2). The 
phase and bandpass calibrator was J0359+509, with boot- 
strapped flux density of 9.5 Jy. We calibrated and imaged 
the data using the standard procedures of the Astronomi- 
cal Image Processing System (AIPS) of the National Radio 
Astronomy Observatory (NRAO). Continuum emission was 
subtracted using task UVLIN. Mapping and deconvolution 
was carried out with task IMAGR, using a natural weight- 
ing of the visibilities. To improve the signal-to-noise ratio for 
extended emission we convolved the resulting images with a 
Gaussian beam, to obtain a final synthesised beam size of 
8". 

We have also reprocessed archival radio continuum data 
of this region at 2 and 3.6 cm. These observations correspond 
to VLA projects AC240 (1989 March 19, B configuration), 
AT121 (1991 May 25, D configuration), and AT122 (1991 
January 30, CnD configuration), and were published in Tor- 
refies et al. (1992; projects AT121 and AT122) and Kurtz et 
al. (1994; project AC240). The description of those observa- 
tions can be found in these papers. Further self-calibration 
was possible for the data observed in 1991. We combined to- 
gether the visibility data at each wavelength before imaging. 
With this combination we can properly sample the emission 
at both, larger (as traced by the Torrefies et al. data), and 
smaller scale (given by the higher-resolution Kurtz et al. 
data). Imaging of the combined data was carried out with a 
robust weighting of the visibilities, with robust parameters 
and 5 for the data at 3.6 and 2 cm, respectively. With these 
parameters, the final data have a similar size of their syn- 
thesised beams (1.01" x 0.84" at 3.6 cm, and 1.18" x 1.11" 
at 2 cm). 

We have also retrieved publicly available submillimetre 
data at 450 and 850 /im, taken with SCUBA at the James 
Clerk Maxwell Telescope (JCMtQ on 2003 November 6, 

^ The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 

^ The James Clerk Maxwell Telescope is operated by the Joint 
Astronomy Centre on behalf of the Science and Technology Fa- 
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Figure 2. Spectra of the NH3(1,1) and (2,2) lines at the position 
of the peaks of the clumps shown in Fig. [T] 

corresponding to proj ect M03BU13. The data at 850 pm 
have been presented inlCurran fc Chrvsostomoul (|2007 ^') and 
iMatthews et al.l l|2009l ). but in this paper we are using for 
our analysis the original data (~ 14" angular resolution), 
without applying any spatial convolution (final angular res- 
olution of 20" in Matthews et al. 2009). The data at 450 /xm 
are presented here for first time. The angular resolution of 
the JCMT at 450 /xm is ~ 8". 



3 RESULTS 

3.1 Ammonia observations 

Fig. [1] shows the integrated intensity maps of the main hy- 
perfine component of NH3(1,1) and (2,2) lines. No signifi- 
cant emission is present at the central cluster of stars. Two 
clumps are clearly present in both maps, to the east and 
south of the centre (labelled as CI and C3 in the figures). 
Their morphology appears elongated, specially in the (1,1) 
map. The spectra toward the peaks of those clumps are 
shown in Fig. [S] The main and inner satellite components 
of the electric quadrupole hyperfine structure of the (1,1) 
transition are evident. A weaker clump (labelled C2) is also 
visible in the (1,1) map, west from CI. In tables [T] and [2] 
we list the observed and the derived physical parameters, 
respectively, of those three clumps. Line intensities, central 
positions and widths in Table [T] were obtained with Gaus- 
sian fits to the emission. Physical para meters in Table [2] wer e 
obt ained following the for mulation of IHo fc TownesI l|l983h . 
and iMangum etahl (Il992l ') . 

Clump CI (the most intense one in NH3) appears to 
have the highest column density. Moreover, its lines are 
wider than the rest. There is an infrared source, detected 
in Spitzer IRAC bands (sll, Kumar Dewangan & Anan- 
darao 2010), close to its maximum, at the exact position 
of the NH3(2,2) peak. These characteristics may imply that 

cilities Council of the United Kingdom, the Netherlands Organi- 
sation for Scientific Research, and the National Research Council 
of Canada. 
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Table 1. Observed properties of NH3(1,1) clumps. 



Clump 


R.A.I 


Dec 1 


1(1, 1; m) 2 


/(l,l;a) 3 






/(2, 2; m) ^ 


V2 ' 


AV2 » 




(.J2000) 


(J2000) 


(m.Ty boam-1) 


(m.Iy boam-1) 


(km 


(km 3-1) 


(m.Ty boam-1) 


(km s-1) 


(km a-1) 


CI 


03 07 27.9 


58 30 34 


53.3 ±1.1 


18.3 ± 1.1 


-41.11±0.02 


2.21 ± 0.05 


23.4 ± 0.6 


-41.04±0.03 


1.96±0.07 


C2 


03 07 25.7 


58 30 42 


28.8 ± 1.0 


12.2 ± 1.1 


-40.00±0.02 


1.29 ± 0.05 


17.6 ± 1.0 


-40.19±0.03 


0.85±0.06 


C3 


03 07 21.4 


58 30 06 


36.7 ± 1.1 


13.1 ± 1.7 


-40.43±0.02 


1.34 ± 0.05 


18.7 ± 1.1 


-40.42±0.03 


1.05±0.07 



1. Coordinates of the positions of the local peak in the integrated intensity map of NH3(1.1) (see Fig. 0. Units of right ascension are hours, minutes, and 
seconds; units of declination are degrees, arcminutcs, and arcseconds. 

2. Intensity of the main hyperfine component of the NH3(1,1) line, at that position. 

3. Intensity of the inner satellite component of NH3(1,1). 

4. LSR velocity of the NH3(1,1) line. 

5. Full width at half maximum (FWHM) of the main hyperfine component of the NH3(1,1) line. It was derived from a simple Gaussian fit, without considering 
the magnetic hyperfine structure within the lines. 

6. Intensity of the main hyperfine component of the NH3(2,2} line. 

7. LSR velocity of the NH3(2,2) line. 

8. FWHM of the main hyperfine component of the NH3(2,2) line. 
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Figure 3. Channel maps of the main hyperfine component of the NH3(1,1) (left) and NH3(2,2) (right) lines. The lowest contour level 
and the increment step are 3 times 3.3 mjy, the rms of the maps (synthesized beam= 8", shown at the bottom left-hand corner of the 
channel with velocity —43.1 km s~-^). 



CI harbours a protostar inside. No infrared source appears 
clearly associated to the peak of clumps C2 and C3. 

Channel maps are shown in Fig.[3]for the (1,1) and (2,2) 
lines. We do not see any obvious velocity pattern within the 
region, although there is a hint of redshifted emission to the 
north (see panel at —38.8 km s~^). This would suggest a 
trend (redshifted to the north, and blue shifted to the south ) 
similar to that found in the CO outflow (jComez et al.|[l99^ . 
although the low signal to noise ratio of the emission, apart 
from clumps CI and C2 precludes any firm conclusion based 
on the velocity pattern. 

We note that single-dish observations (obtained with 
the Effelsberg antenna, angular resolution ~ 40"; Wu et al. 
2006) show ammonia emission associated with the central 
cluster, specially in the (2,2) transition. Since we do not 
detect this in ou r observat i ons, i t is possible that the emis- 
sion detected bv lWu et al.l (|2006l ) is extended and relatively 
uniform, so that either it is below the sensitivity limit when 
observed at 8" resolution, or a fraction of it is missed because 
of the lack of short interferometer spacings. To estimate the 
amount of emission that our interferometric observations are 
missing, we have convolved our integrated intensity maps 



Table 2. Physical parameters of NH3(1,1) clumps. 



Clump 








size * 


M 5 






(K) 


(10^2 cm-2) 


(arcsec) 


(M©) 


CI 


0.6 


18 


1.4 


21 X 12 


5 


C2 


1.4 


20 


0.5 


11 X 8 


0.7 


C3 


0.8 


19 


0.5 


40 X 13 


4 



1. Optical depth of the main hyperfine component of the NH3(1,1) 
line. 

2. Rotational temperature. 

3. Column density of hydrogen, assuming local thermodynamic 
equilibrium, with kinetic and excitation temperature of NH3 lines 
equal to Tr, and an abundance of NH3 relative to hydrogen of 
10-8. 

4. FWHM size of the clump. 

5. Estimated mass of the clump. 

with an elliptical Gaussian, to obtain a final angular res- 
olution of 42.2", corresponding to that of Effelsberg at the 
ammonia frequency. We then compared the integrated inten- 
sities tabulated by Wu et al. (2006) with the values in our 
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convolved map. At the position given by Wu et al. (2006), 
the VLA is recovering ~ 12% and 33% of the NH3(1,1) and 
(2,2) emission, respectively. Therefore, we are missing a sig- 
nificant fraction of the NH3 emission, which would compli- 
cate any conclusion based on large-scale structures. How- 
ever, at smaller scales, we are confident that the individual 
ammonia clumps seen in Fig. 1 are real structures, given 
their good correspondence with submillimeter emission (see 
next section). To get a complete picture of the overall dis- 
tribution of ammonia map, we should combine both inter- 
ferometric and single-dish data. Unfortunately, the available 
single-dish data are not enough to properly carry out such a 
combination, since the Wu et al. (2006) maps are not fully 
sampled (their observed points have full-beam separation). 

3.2 Submillimeter continuum emission 

The SCUBA maps a t 85 /xm presented in 
ICurran fc Chrvsostomoul |2003) and [Matthews et al.l 
1 20091 ) showed extended, elongated emission oriented in 
the northwest-southeast direction. However, in the maps 
at 450 /xm, (Fig. [4]), we can clearly distinguish two local 
maxima. One of them is located at the position of the 
central cluster, and the other one coincides with CI (with 
its maximum located at R.A.(J2000) = 03''07'"28.1^ 
Dec(J2000) = 58°30'35")- We also note that the central 
submillimeter emission extends to the east, toward C2. 
There is also a weaker finger of emission elongated toward 
the south that peaks at the position of the IR source s4 
(see Fig. and extended, very faint emission coincident 
with C3. Apart from the central peak, the submillimeter 
emission tends to follow the distribution of NH3 . This gives 
us confidence that the ammonia clumps seen in Fig. [T] are 
real entities, and not artifacts generated by the lack of short 
spacings in the interferometer. 

From the map at 450 fim, we can estimate the mass 
of the observed submillimeter clumps, using the formula 
M = Svd? /{KvBi,{Td)), where Si, is the fiux density, d is 
the distance to the source (2 kpc), Hi, is the mass absorption 
coefficient, and B,^{Td) is the Planck function for the dust 
temperature, Td- The value of is highly uncertain, and is 
very sensitive to the particular dust properties o f the object. 
We ha ve adopted the prescription suggested by iHildebrandl 
(|l983l ), K„ = 0.1(250Mm/A)'' cm g , which assumes a gas 
to dust ratio of 100. To estimate the dust emissivity index 
{(3) we convolved the 450 /im map with an elliptical gaus- 
sian, to match the angular resolution of the one at 850 fim. 
After this convolution, the ratio between those maps is ~ 9, 
which yields /3 ~ 2.2, a value within th e range found in other 
high-mass stars l|Molinari et al.l [ 20001 ). With these parame- 
ters, fiux densities of 6.9 and 5.6 Jy represent masses of 44 
and 39 Mq for the central submillimeter peak and the one 
associated with the ammonia clump CI, respectively. In the 
latter case, the mass obtained from the submillimeter data 
is larger than the one from ammonia (5 M0), but consistent 
with our estimate that the ammonia data may be recovering 
only ~ 12% of the total fiux. 

3.3 Centimeter continuum emission 

Maps of the continuum at 3. 6 and 2 cm are shown in 
Fig. [5] As reported previously l|Wvnn- Williams et al.|[T98ll : 
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Table 3. Centimeter radio continuum sources. 



Source 


5(3.6cm) 1 


5(2cm) 2 


a(3.6 - 2) ^ 




(mJy) 


(mJy) 




W 


22.1 ±0.7 


28.3 ± 2.2 


0.4 ±0.2 


S 


1.5 ± 0.4 


< 2 


< 0.5 



1. Flux density at 3.6 cm. 

2. Flux density at 2 cm. 

3. Spectral index between 3.6 and 2 cm, defined as S,j (x v°' . 

iTorrelles et al.|[l99^ : iKurtz et al.lll994h . there are two main 
radio continuum sources, as seen in Fig. [5] associated with 
sources AFGL 437W and S. Parameters of the continuum 
sources are shown in Table [3] However, our map at 2 cm 
shows extended, elongated emission in AFGL 437W, in the 
NE-SW direction (P. A. ~ 60°). This elongated emission was 
not clearly seen in the individual datasets that we combined 
together to obtain this map. It is possible that the elongated 
structure only shows up with good enough coverage of the 
uv-plane, which would also explain why it is not seen in the 
3.6 cm data. If this interpretation is correct, further obser- 
vations at 3.6 cm with a better uv-coverage could detect the 
elongated structure. 

The elongated radio continuum emission suggests the 
presence of a jet, and it is reminiscent of other jets traced 
by radio continuum in high-mass star-forming regions, like 
Cepheus A or W75N (Torrefies et al. 1996, 1997). The de- 
rived spectral index of this source (0.4 ± 0.2) is similar to 
that expected from constant-velocity winds (~ 0.6). Follow- 
ing the formulation by Rcvnolds (1986), the spectral index 
would indicate the presence of a confined jet (i.e, narrower 
than a purely biconical shape). This formulation assumes a 
power-law dependence of width (w) with distance from the 
center (r), of the form w cx r^. In our case, the derived index 
is e ~ 0.8. 

The mass loss rate of the ionized jet can be estimated 
from eq. 19 of Reynolds (1986). We assumed a wind velocity 
of 500 km s~^, a turnover frequency for the radio continuum 
emission of 10 GHz, a temperature of the ionized gas of 10* 
K, a jet opening angle of ~ 30°, and an inclination angle 
with respect to the line of sight of 45°. We have chosen 
a turnover frequency between the two observed ones, since 
the derived spectral index (0.4) is partially optically thick. 
However, the value of the mass loss rate depends only weakly 
on the turnover frequency (to the power of —0.15), and the 
possible error introduced by this value is small. Under these 
assumptions, we obtain a mass-loss rate of ~ 5 x 10~® M© 
yr~^. The assumed velocity of 500 km s^^ would imply a 
momentum rate of ~ 2.5 x 10""^ Mq km s~^ yr"'^ for this jet, 
which could be enough to drive the large scale CO outflow 
in the region, given that the derived momentum rate of jet 
and outflow are of the same order (Gomez et al. 1992). 

We also note t hat mid-infrared images of this source 
(Ide Wit et al .l l2009l ) show extended emission slightly elon- 
gated along the same direction. 

In the 2 cm map, there seems to be also some radio con- 
tinuum emission ~ 3.7" northeast of AFGL 437S, and an ex- 
tension to the southeast of AFGL 437W. We have searched 
in data archive for possible infrared counterparts of those 
emissions, but none was found. At this point, we cannot 
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Figure 4. (Left) Map of the emission at 450 fim. (Right) An overlaid of 450 /im emission (contours) with integrated emission (greyscale) 
of the main component of the NH3(1,1) line. Contours are -5, 5, 10, 15, 20, 25, 30, 35, 40, 45 and 50 times 0.152 Jy beam~^, the rms of 
the map (beam~ 8", shown at the bottom left-hand corner of the right map). 
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Figure 5. (left) Radiocontinuum emission at 3.6 cm in the central cluster of AFGL 437. Contour levels are -4, 4, 6, 10, 15, 20, 25, 30, 
40, 50 and 60 times 0.042 mjy beam~i, the rms of the map (beam= 1.01" x 0.84", shown at the bottom left-hand corner of the map). 
(Right) Radiocontinuum emission at 2 cm. Contours are -4, 4, 5, 6, 8, 10, 12, 15, 20, 25, 30, 35 and 40 times 0.1034 mJy beam "i, the 
rms of the map (synthesized beam= 1.18" x 1.11"). 



ascertain whether these represent two real sources, or they 
are instrumental artifacts. If they are indeed real sources, 
the lack infrared counterpart for the weak 2 cm continuum 
emission may suggest the presence of young, deeply embed- 
ded sources. Some weak continuum emission may also be 
associated with AFGL 437N. 



4 DISCUSSION 

Gomez ct al. (1992) suggested an interpretation for the 
molecular outflow in terms of motions of gas along the wall 
of a cavity opened in the molecular cloud by the massive 
stars in the cluster. With this anisotropic distribution of 
molecular gas, even if the stellar winds are isotropic, the 
outflow could appear as bipolar, and of low coUimation, to 
the observer. 
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Our NH3 observations were specifically designed to test 
whether this "alternative" interpretation for the outflow is 
valid, or this is a more "classical" bipolar outflow, where 
there is intrinsic coUimation in the mass-loss process, prob- 
ably driven by source WK 34. There are several key char- 
acteristics to be investigated with these observations, such 
as the velocity pattern, the morphology, and the location of 
NH3. 



distinctive nature as a submillimeter source shows that sll 
is a YSO. Nevertheless, source sll is clearly offset from the 
center of the outflow, and it is hard to imagine that it could 
be its main driving source. However, we note that both, the 
redshift ed and blue shifted CO lobes show extensions to the 
east (iGomez et al.|[T992l . We suggest that these extensions 
could in fact be part of an independent outflow driven by 
source sll, while the bulk of the CO outflow would be driven 
by the sources in the central cluster. 



4.1 The velocity pattern of NH3 

If the interpretation given by iGomez et~aLl l|l992l) for the 
outflow, as motions along the walls of a cavity is correct, 
this drag of material should be seen in all molecular tracers. 
In particular, the spatial distribution of blue- and redshifted 
velocities of the dense gas traced by NH3 should be similar 
to the one seen in CO. 

Our results (Fig. (2]) are not conclusive to reject the 
Gomez et al.'s hypothesis. There is some hint of more red- 
shifted gas to the north, and more blueshifted one to the 
south, but the low signal-to-noise ratio of the emission, spe- 
cially that to the north does not allow flrm conclusions. 
On the other hand, the most intense ammonia emission 
(CI) seems to be related to a particular infrared source 
(sll), which also emits in the submillimeter. The ammonia 
seems to trace individual clumps, well defined in velocity, 
rather than a continuous distribution with a velocity gra- 
dient, which one co uld expect from the model proposed by 
iGomez et all (|l992l ). 



4.2 The morphology and location of NH3 

Classical examples of bipolar outfiows show interstellar 
toroids of dense gas, perpendicular to molecular outflows 
and/or jets (e.g., Torrelles et al. 1983; Wiseman et al. 2001). 
In the past, these interstellar toroids have been proposed as 
the coUimating agents of these outflows, although the ac- 
tual collimators seem to be much smaller circumstellar disks. 
However, given the usual relationship outflow-interstellar 
toroid, the presence of a dense structure (traced with NH3) 
perpendicular to the outflow, and located close to is cen- 
ter, would have favored a more "classical" interpretation 
of the outflow. Moreover, ammonia emission has been used 
as a tool to identify the excitation source of bipolar out- 
flows. Excitation source of outflows usually coincides wit h 
the maximum emission of ammonia (|Anglada et all 1 19891 ). 
This is usually confirmed with the presence of local enhance- 
ments of temperature (with a higher ratio of the emission of 
ammonia (2,2) to (1,1) transitions), and turbulence (wider 
ammonia lines). 

The overall distribution of NH3 does not show any pref- 
erential orientation with respect to the molecular outfiow. 
Only clump 3 shows an E-W orientation, roughly perpen- 
dicular to the outflow. However, it is located to the south, 
on the blueshifted lobe of the CO emission, rather than to- 
wards the central cluster, and no infrared source is located 
near its maximum. 

Using the maximum of ammonia (clump CI) as a cri- 
terium to search for the powering source of the outflow, 
source sll w ould be a good candidate. Its spectra l energy 
distribution jKumar Dewangan k, AnandaraolbOlOl ) and its 



4.3 Sources of mass loss in the region 

There are several sources in AFGL 437 that could be under- 
going mass-loss simultaneously. Sour ce WK 34 is associated 
with an infrared polariz ed nebula (|Weintraub fc KastneJ 
ll996l : lMeakin et akllioosl ). and has been suggested to be the 
dominant exciting source of the outflow. The inf rared neb- 
ula and its association with a wa ter maser (iTorrelles et al.l 
ll992l : IWeintraub fc Kastnerlll996l 'l clearly indicate that this 
source is an active source of mass-loss. Source AFGL 437 S 
also shows elongated infrared emission (| Alvarez et al.|[2004l '). 
which is also suggestive of a star undergoing mass-loss, al- 
though this extension is weak, and close to the resolution 
limit of the maps. 

Our radio continuum map at 2 cm (Fig. [S]) shows, for 
the first time, the presence of a coUimated jet associated to 
source AFGL 437W. The previous detection of both radio 
continuum and water maser emission (Torrelles et al. 1992) 
already signaled this source as a young object, but now we 
see that its mass loss is highly coUimated. 

Therefore we have strong evidence that at least WK 34 
and AFGL 437W (and possibly AFGL 437S) are undergoing 
co Uimated mass loss. This seems not to support the model 
of lGomez et"aLl (|l992l ). which suggested that the mass loss 
could be isotropic (and the observed CO bipolarity would be 
a projection effect). We now favor that the observed CO out- 
flow is the superposition of several bipolar outflows, which 
would explain its low degree of coUimation. Sources WK 34, 
AFGL 437W (and possibly AFGL 437S) could be respon- 
sible for the bulk of high-velocity CO emission close to the 
center. We also propose that source sll is a young embedded 
object that could drive an additional outflow, traced by the 
eastern ext ensions of the CO lobes. Weintraub fc Kastiierl 
(1996) andU umar Dewangan fc Anandarad l|2010l ) noticed 
that the infrared emission near WK 34 is oriented in the N-S 
direction next to the source, but it bends to the nort heast 
away from it. iKumar Dewangan fc Anandarad (|2010l ) sug- 
gested that this bending may be due to the interaction of 
the outflow from WK 34 with mass loss from source AFGL 
437N, but we suggest that it could reflect the interactions 
between the outflows from WK 34 and AFGL 437W. 

This superposition of several individual outflows, giving 
rise to complex morphologies when observing at low angular 
resolution, is also found in other high-mass star forming re- 
gions (e.g., Beuther et al. 2002, 2006; Ginsburg et al. 2009). 

An alternative scenario would be that only one source 
is signiflcantly driving the (low-coUimation) outflow, even 
when present evidence is of highly coUimated jets, if both a 
coUimated jet and a low-coUimation wind are driven simul- 
taneously by the same source. Evidence for simultaneous 
presence of high- and low-coUimation mass loss has been 
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found for the first time in a high- mass youns star in w ater 
maser observations of Cepheus A (jTorrelles et al.ir201lh . 

To test these scenarios, we propose that the molecular 
outflow should be observed with millimeter/submillimeter 
interferometers, at scales of a few arcseconds, to be able to 
discriminate the possible individual outflows driven by these 
sources. Being able to trace the molecular outflow close to its 
driving source(s) would provide direct evidence of whether 
there is a superposition of outflows or if one source is re- 
sponsible for it. 



5 CONCLUSIONS 

We presented new interferometric NH3(1,1) and (2,2) obser- 
vations of the region of high-mass star formation AFGL 437, 
with the aim to understand the nature of the low-collimation 
bipolar molecular outflow previously detected in CO (Gomez 
et al. 1992). These observations were complemented with 
archive data of radio continuum emission at 2 cm, 3.6 cm, 
and 450 yLxn. Our main conclusions are as follow: 

(i) The ammonia emission in our interferometric maps is 
located mainly in three clumps at the south and east of the 
central cluster. 

(ii) The analysis of the velocity pattern of the ammo- 
nia transitions in the whole region is not conclusive and, in 
principle, it does not allow us to accept or reject the model 
proposed by Gomez et al. (1992). 

(iii) One of the clumps seen in NH3(1,1) coincides with 
the ma:ximum of NH3(2,2) emission, and with a local maxi- 
mum of emission at 450 fj,m. A near infrared source (sll) is 
also found at that position. We identify this source as pos- 
sible active young stellar object, and it might be a source of 
mass-loss in the AFGL 437 region. 

(iv) The radiocontinuum map at 2 cm shows extended, 
elongated emission in the source AFGL 437W, not detected 
in previous studies. This elongated morphology and its spec- 
tral index (~ 0.4) suggest the presence of a jet in that source. 

(v) The coUimated mass loss observed in sources WK 34 
and AF GL 437W sour c es do es not support the model pro- 
posed by I Gomez et al.l (|1992D . which suggested an isotropic 
mass loss in the region. 

(vi) We suggest that several bipolar outflows may exist 
in the AFGL 437 region. The CO outflow would be the 
superposition of them, what would explain its low degree 
of coUimation. An alternative scenario is that the CO out- 
flow could be produced by a main source driving both less 
and highly coUimated winds. Observations with millime- 
ter/submillimeter interferometers to scales of several arc sec- 
onds may identify the individual outflows. 
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